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Abstract 

The statistical properties of overlap sums of groups of four eigenfunctions of the Anderson model for localiza- 
tion as well as combinations of four eigenenergies are computed. Some of the distributions are found to be scaling 
functions, as expected from the scaling theory for localization. These enable to compute the distributions in 
regimes that are otherwise beyond the computational resources. These distributions are of great importance for 
the exploration of the Nonlinear Schrodinger Equation (NLSE) in a random potential since in some explorations 
the terms we study are considered as noise and the present work describes its statistical properties. 



1 Introduction 



In the present work the statistics of the overlap sums and combinations of the eigenenergies of the one dimensional 
Anderson model [HE] are calculated. These quantities naturally arise in the exploration of the Nonlinear Schrodinger 
Equation in a random potential [3j [4J [5j El [JJ [8] . The overlap sums defined by ^ measure the overlap of four 
eigenfunctions of the Anderson model. The combinations of the eigenenergies dominate the phase of the nonlinear 
terms. As a result of the nonlinearity such terms that affect the dynamics, play a crucial role in the dynamics of 
the nonlinear model. Sometimes these terms are considered as effective noise [HI [5]. Their statistical properties 
are crucial for the effective noise theories. The statistics presented here may be also of interest for mathematicians 
exploring the Anderson model and related random models. The Nonlinear Schrodinger Equation (NLSE) in a 
random potential takes the form [3] [31 El [7J |SJ E] 



id t ^ = H^ + (3\^\ 2 ^, 



(1) 



where Hq is the linear part with a disordered potential, which on a lattice takes the form of 



H ip(x) = - (ip(x + 1) + ip{x - 1)) + s(x)i/j(x). (2) 

In this work it is assumed that e (x) are identical independent random variables (i.i.d) uniformly distributed in the 
interval of [— • This is the famous Anderson model [TJ[5]. Here we study the model in one dimension where 
all the states are exponentially localized (TUI ITTj 

The NLSE was derived for a variety of physical systems under some approximations. It was derived in classical 
optics where ip is the electric field by expanding the index of refraction in powers of the electric field keeping only the 
leading nonlinear term |12[ PTSj fT4"] . For Bose-Einstein Condensates (BEC) , the NLSE is a mean field approximation 
where the term proportional to the density /JIV'I 2 approximates the interaction between the atoms. In this field the 
NLSE is known as the Gross-Pit aevskii Equation (GPE) [Tg| [T7| [T51 [TO] . 

A natural question in this research is whether a wave packet that is initially localized in space will indefinitely 
spread for dynamics controlled by [3J HI |5J [El IS] ■ 

It is convenient to expand the wavefunction ip in the linear problem cigenfunctions u n 

^(x,t)=Y / Cn(t)u n (x)e' lE " t (3) 

n 

where E n is the eigenenergy of the n-th eigenstate. The width of the energy spectrum is A = 4 + W with 
E n 6 [—2 — 2 + ^1 and u n typically falling off exponentially [101 [TT] 

e -|x n -x|/£„ 

u n (x) ft = ip(x) (4) 



where <p{x) is a real random function of order unity, £„ is the localization length for the n-th state and x n is its 
localization center. The expansion coefficients c n (t) satisfy pUH] 



id tCmi {t)=P E ^ 2 ' m3 ' m4ex P( i -* $ mr m3 ' m4 )C 2 Cm 3 C m4 (5) 
m2,m 3 ,m 3 

where 



is the overlap sum and 



Cr 3 '™ 4 =J2u mi (x)u m2 (x)u m3 (x)u mi (x). (6) 

x 

fom 2 ,m 3 ,m 4 = ]?+E~E~E (7) 

is the total phase. 

The main purpose of this paper it to explore numerically the statistical properties of the overlap sum yw>> m 3> m * 
and the total phase $™2> m 3> m 4 as well as of related quantities. We believe that understanding the qualitative 
properties of V™f' m3 ' mi and $™2,m3,ni4 w ju ne jp ^ build a toy model that may shed light on the spreading 
mechanism. Of particular interest will be to find universal distributions in terms of scaling variables as explained 
in Sec. 2. These should be relevant for the entire energy regime. The distributions are found numerically and are a 
generalization of the distribution of V^ ' ' that was found analytically by Fyodorov and Mirlin in a narrow energy 
range j2H|- In Sec. 2 we explore statistical properties of V^' m3 ' m * for weak disorder and in Sec. 3 we will explore 
statistical properties of <j>™2,m3,m4 f or wea k disorder. For weak disorder it is expected from the scaling theory that 
all statistical properties are determined by the localization length, while for the strong disorder this does not hold 
[21J. Some results for strong disorder are presented in Sec. 4. The results are summarized and discussed in Sec. 5. 
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Statistical properties of VJ^ 2,m3,m4 for weak disorder 



mi 

We explore the values of y^ 2 '" 13 - 7 ™ 4 j n the regime of weak disorder (which corresponds to relatively long localiza- 
tion length £). Most of the explorations are numerical. The lattice size is fixed at N = 500. For each realization of 
the e (x), we computed the eigenfunctions u n and ordered them in space by the center of norm coordinate, defined 
by x n = J2 X x ' u n ( x ) ■ We choose u mi with mi = to be the eigenfunction centered in the lattice. We studied only 
the following quantities V^ ' 0,0 , V ' ' , V 0,0,1 ,V 0,1,2 and V 12,3 representative of values where y™2,m 3 ,m4, are i ar g e 
Combinations with mi > £ (taking m\ — 0)have negligible values because the overlap sum is a sum of exponentially 
decaying functions in space of the form of Q. We calculated these values for Nr — 2- 10 4 realizations, and repeated 
this calculation for 7 disorder strengths in the weak disorder regime 1 < W < 2 where the maximal localization 
length £ takes the values of 25 < £ < 103. In the regime of weak disorder the maximal localization length is [22] 

We computed the distributions of the j/™ 2 *" 13 ' 1 ™ 4 as follows. We calculated 2 ■ 10 4 values of V^ m3 ' mt , one for 
each realization. We know that jv^ 2 '" 13 '" 1 - 4 1 must satisfy < |]/™2.m 3: m 4 | < -y ^^^^ t ne eigenfunctions are all 
normalized ^ x it 2 (x) = 1. We made a histogram of the values jy^ 2 '™ 3 '™ 4 1 j n number of bins Nbina = 500 in the 
interval [0, 1], the resulting bin size is Sx = 0.002. In order to get the distribution we normalized the values of the 
histogram, dividing them by the number of realizations, Nr. 

In the calculation of the statistical properties of the y^ 2 '™ 3 '" 14 we distinguish different groups according to the 
number of different indices m,\ 

2.1 The case mi = ni2 = = = 

In the case where all indices are equal we have chosen them to be zero. In this case 

<°'° = E M o(-)^o. (9) 

x 

It is just the inverse participation ratio. Its distribution was calculated analytically by Fyodorov and Mirlin |20| 
and was found to satisfy scaling, that is if P (Vo, £) is the probability density of Vq and the localization length is £ 
then, if one defines a scaling variable 

yo = v £ (10) 

its probability density is 

P(yo) = ^P(V ,Z). (11) 

In [20J this scaling was found to hold in a narrow range of energy. In the present work we demonstrate numerically 
that it is an excellent approximation also when the maximal localization length Q is used. The scaling function is 
different from the one of |20| . 

First we verify that the average of Vo satisfies 

(Vo) = j (12) 
where C is a constant independant of £, and £ is given by (pi, as maybe expected from the scaling relation (in 



3 



agreement with [5]). This is clear from Fig. [T] and it is found that C — 1.296 . . .. 

The probability density function (PDF) as a function of Vq is presented in Fig. [2Jl . A typical function fitted to 
the numerical data is shown in Fig. [2J3 for W = 1 and it takes the form 

P(l/ ) -c ie - (c2+C3 - ln(y ° ))2 (13) 



with ci = 127.7 . . ., C2 = 9.097 . . . and C3 = 1.865 .... We found that the scaling ( 10 1 and (111 holds for all weak 
disorder strengths studied as shown in Fig j2jj. The resulting scaling function is 

P (y Q ) = a ie - (Q2+a3 ln(yo))2 . (14) 

with ai = 1.21 a 2 = 0.539 . . . and a 3 = 1.71 . . .. 

What is the reason for the scaling? From ^ and ^ it is clear that the magnitude of each of the u mi is of 
order ^= while the number of terms in the sum that contribute substantially is of order £ . Therefore Vo, although 

random, it is typically proportional to 7^-. Note that all the contribution to the sum Q 6 p are positive. 



If the calculation is confined to a narrow energy, £ is practically constant and P (j/o) is the function found in 
|20| . In the case we study the energy of the site mi = (middle of the lattice) varies as the realizations change and 
an effective average over the realizations is preformed. Since the density of states (see Fig. |6| and the localization 
length as a function of energy are flat at the center of the band, where the localization length is maximal and takes 
the value close to ^ , terms with this value of the localization length dominate the overlap sum (J6j> . It is worthwhile 



to note that the scaling function ( 14 1 we found is different from the one found in [20 . It is practically the average 



of the function found in [20 over energy. 

Now we consider the cases where the vtii take two different values say V ' ' of Vq ' ' 1 . 

2.2 Distribution of V ° 1 1 

Also here an argument similar to the one presented in the previous section holds, but the localization lengths of 
the two wave functions involved are different, the overlap sum is of the order ^ ^_c x , therefore Vq ' 1 ' 1 behaves as |. 

In order to investigate the distribution of Vq ' 1 ' 1 , denoted here by V\, which consists of many near zero values, we 
generated the histogram of In (Vi) rather than V\ . In Fig. |3]we present the distribution of P (In (Vi)) as a function 
of In (Vi). The best fit for the scaling function, in terms of the scaling variable 

Vi = ViZ (15) 



is shown there as well. As expected (Vi) satisfies a relation similar to (12) but with C = 0.429. . . (in agreement 
with 0). 

2.3 Distribution of V 1 

Let us denote Vq ' ' 1 = V^. From the definition ^ it is clear that (V2) = 0. Therefore to estimate the typical 
value of V2 we study ( V 2 2 ) . It can be estimated by 
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Figure 1: a log-log scale of z = In ((Vb)) as a function of x = In (£) . The numerical results are represented by circles 
while the line is the best fit. The fit is z — —a$x + bo for ao = 0.95 . . . and 60 = 0.26 . . .. 

It is of order r2 ? s ■ Therefore it is reasonable that (V 9 2 ) ~ Indeed one finds 

yfty?) = cc 1 -* (17) 

with C = 0.566 . . . independent of £. This motivates us to introduce the scaling variable 

y2 = v 2 e- 5 . (is) 

In Fig. |4]we show the distribution of P (In (V^ 2 )) as a function of In (V2) in terms of the scaling variable y 2 - 
2.4 Distribution of v^' m3 ' m4 when 3 or 4 different m; are involved 

In this case we could not find any simple scaling relation. The averages are found to be exponential in £, as one 
can see from Fig. [5] 

3 Statistical properties of $^i> m 2> TO 3 f or wea k disorder 

In this section we explore the statistical properties of $ defined in Sec. 1 and the distribution of the eigenenergies 
E n . For the weak disorder regime we fixed the lattice size N = 500 and computed the eigenenergies for Nr = 10 3 
realizations. We repeated this calculation for 7 different values of the disorder strength ,1 < W < 2 which correspond 
to 25 < £ < 103. 
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Figure 2: The PDF for various values of weak disorder strength, (a) The PDF as a function of Vb for: blue circles 
W = 1 (£ 103), green squares W = l\ (f « 75) , red crosses W = 1§ (f « 58), turquoise dots VF = l| (f « 46), 
purple pluses W = if w 37), green stars = if (f w 30), black rhombus PF = 2 w 25). (b) The PDF for 
W = 1, the results of the simulation are presented by blue circles, while the solid line is given by (131. (c) P (yo) 



as a function of y . The data collapse indicates the scaling of the PDF with the localization length £. The colored 
lines correspond to various values of the disorder strength. The black thick line is the best fit for the log normal 



distribution (14) 
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Figure 3: The PDF of ln(Vi) for various values of weak disorder strength presented in Fig [2] (a) The PDF as a 
function of In (Vi) (b) The scaled PDF as a function of In (yi) the black thick line represents the best fit. The data 
collapse indicates the scaling of the PDF with the localization length £. The colored lines correspond to various 
strengths of disorder. 
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Figure 4: As in Fig. |] the PDF of In (V 2 2 ) for various values of weak disorder strength presented there, (a) The 
PDF as a function of In (V^ 2 ) . (b) The scaled PDF as a function of In (y 2 )- The black thick line is the best fit. The 
number of bins here is 200. 

3.1 Distribution of E n 

The distribution of the eigenenergies for the weak disorder regime, as plotted in Fig. [6] is symmetric around 
E = and characterized by convex function in the middle and sharply decaying function at the boundaries. 

3.2 Distribution of $+ = E n + E m 

We calculated the distribution of the sums of two eigenenergies obtained for the same realization. The motivation 
for calculation of these sums is from the terms where m\ = m,2 = and 7713, are arbitrary in j5|. In Fig. [7] we 
plot distributions of $ + with various disorder strengths. Note the maximal value of the distribution decreases with 
W. We found the following relation between the maximal value of the distribution which we will denote <& m and £, 

$ m (0 = 0.13C 015 (19) 

3.3 Distribution of $ " = E n - E m 

We calculated the distribution of the differences between two eigenenergies obtained for the same realization. 
Despite the symmetric nature of the distribution of Ei the value of $~ differs from $ + because of level repulsion 
[2"5l 121] . Therefore we anticipate a relatively large peak at $~ = 0, as shown if Fig. [8j A comparison between i> _ 
and $ + is presented in Fig. [9] One can see that, the distributions differ substantially only near <E> - = and $ + = 
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Figure 6: The PDF of E n for various weak disorder strength, blue circles W = 1 (£ « 103), green squares W = lg 
(£ ~ 75) , red crosses W = l| (£ ss 58), turquoise dots W = \\ (£ « 46), purple pluses T'F = l| (£ « 37), green 
stars W = If (£ w 30), black rhombus W = 2 (£ w 25). 
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Figure 7: The distribution of <5 + for strengths of disorder as in Fig. fusing the same symbols, 
because of level repulsion. 
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Figure 8: (a) The distribution of <I> - for strengths of disorder as in Fig. [6j Number of bins used is 
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Figure 9: (a) The PDF of $ + (blue squares) and $ (green circles) for W = 1, total number of bins used is 10 2 . 
(b) zoom of (a) around <£>+ = and $~ = 0, total number of bins used is 10 3 . 



3.4 



Distribution of ^,m 3 ,m 4 



$™J : " i3 ' m4 = $ is a combination of 4 eigenenrgies. In order to calculate the distribution of $ one needs to 
compute N A of these combinations. To avoid lengthy computations we are presenting a much smaller number 
of realizations and use a smaller lattice size. In this subsection we will present these distribution for Nr = 10 
realizations on a lattice with size N = 128 for 7 disorder strengths in the range 1 < W < 4, which correspond to 
6.5 < £ < 103. A Gaussian like distribution is found for all values of W, as shown in Fig. [TUJl. The form of the 
distribution is 

P($)=Ae^ (20) 

where A is the normalization constant and a is the width of the gaussian. A fit is presented in Fig. [Top for W — 1 
with the values of A = 0.1335 . . . and a = 4.278 ... in agreement with |J. Next we calculated the width of each 
Gaussian a as a function of £ and found 

<j(0=bi-C b2 +b 3 (21) 



with bi = 5.924 . . ., &2 = 0.865 . . . and 63 = 4.17 . . ., this function is presented in Fig. [TOp . For the case of very 
weak disorder, £ ^> 1, we see the value of a approaches a — > 4.17 . . ., which is in agrement with the value found for 
the distribution plotted in Fig. IlOH where W = 0. 
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Figure 10: The PDF of $ for various values of weak disorder strength, (a) blue circles W = 1 (£ ~ 103), green 
squares W — \\ (£ « 46) , red crosses W — 2 (£ pa 25), turquoise dots W = l\ (£ pa 16.5), purple pluses PF = 3 
(£ 11.4), green stars PF = 3^ (£ pa 8.4), black rhombus W = 4 (£ pa 6.5). Number of bins = 1000. (b) the PDF 
as a function of <& with W — 1, the red dashed line is the fit (20 1, 
(d) The PDF as a function of $ for W = 0, 



the fit [21 
and cr = 4.27. 



(c) u as a function of £, the red dashed line it 
the red dashed line it the fit (20) with values A — 0.1335 
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4 Strong disorder 



In the case of strong disorder one does not expect scaling to work j5T]. Indeed we could not find a scaling 
distribution for Vq and V\ defined in Sec. 2 for the regime of strong disorder. Their averages scale with different 
powers of £ as is clear from Fig. 
Fig " 
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The distribution of the E n exhibits a maximum near E = as one can see from 
while for weak disorder a minimum is found there (compare Fig. [6]to Fig. 12). The distributions of $ + and 
presented in Fig 13 exhibit a linear dependence on the values of <I> + and $~ respectively. The distribution of 
$ is similar to the one found for weak disorder, Fig. [Top fits even better gaussian distribution. 
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Figure 11: (a) z = ln((Vb)) as a function of In (£) with a linear fit z = —0.39 • ln(£) — 0.63. (b) z = (V\) as a 
function of £ with the fit z = —0.01 • £ + 0.13. (c) z = In ((V2)) as a function of £ with a linear fit z = —0.21 • £ — 2 
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Figure 12: The PDF of E n for various strong disorder strength, blue circles W = 6 (£ ps 2.85), green squares VF = 6| 
(f ps 2.3) , red crosses VF = 7§ (£ ps 1.9), turquoise dots VF = 8 (f pa 1.6), purple pluses VF = 8§ (£ ps 1.36), green 
stars VF = 9§ (£ 1.17), black rhombus W = 10 (£ ps 1.02). 
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Figure 13: Distribution function in the regime of strong disorder. Symbols are as in Fig. 
(b) The PDF of 
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(a) The PDF of $ H 



5 Summary and Discussion 



In this work the distribution function of some overlap sums V™ 2 '"* ""' 



. of (j5| were calculated. In some cases of 
weak disorder it was demonstrated that scaling holds. In particular it was shown that for weak disorder (or large 
localization length £) the distribution functions of V^ ' 0,0 , V^ ' 1,1 T ""*"'' 

via one scaling variable, y — ^V^ 0,0,0 ; y\ = £Vn ' ' and y 2 = £ 3 ^ 2 K ,0,1 see dlOb , (|15| and (|18|. The distribution 

rff,1.2 _„j T ^1,2,3 



V D : ' are functions of these variables and of £ 



function of yo was calculated numerically (14). We could not find a scaling function for V ' ' and V ' ' . In 
addition to the fundamental interest, the scaling function can be extremely useful for the case when £ is very large. 
It enables to obtain the distribution of V^' m3 ' m4 for regimes where numerical calculations require a basis of size 



that is beyond the available computer resources. Also the averages and variances of the V. 



rai ,7713 , m 4 



were computed 



for some {m,} . In most cases of weak disorder we analyzed, these scale simply with £ as one could guess from the 
distribution functions of the scaling variable. In some cases the averages are exponential in £ (see Fig. [5|. 

The distribution functions of combinations of energies of the Anderson model were studied as well. A difference 
between $ + = E n + E m and $~ = E n — E m was found near $ + = 0, $ _ = 0. It is a signature of level repulsion. 

E„ l3 ~ E mi was found to be gaussian and the dependence of the 



The distribution of $£g,m 3 ,rr»4 = E mi 



En 



variance on the localization length was computed (21 1 



Some results for the various distributions were calculated in the regime of strong disorder. Results that are 
typically different from ones found for weak disorder were obtained. 

The research can be continued in two directions. The distribution can be calculated numerically for additional 
cases and divided into classes with different characteristic properties. Also better analytical insights are called for. 

This work was motivated by the effective noise theory for Nonlinear Schrodinger Equation (NLSE) in a random 
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potential [HI E [S] and the results of this work are of great importance for this theory. However the results may be 
of interest by themselves. 

We would like to thank Y. Krivolapov for detailed discussions, extremely valuable technical detailed help and for 
extremely critical reading of the menuscript . This work was partly supported by the US-Israel Binational Science 
Foundation (BSF), by the Minerva Center of Nonlinear Physics of Complex Systems, and by the Shlomo Kaplansky 
academic chair. 
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